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Abstract: Brane gas cosmology provides a dynamical decompactification mechanism that 
could account for the number of spacetime dimensions we observe today. In this work we 
discuss this scenario taking into account the full bosonic sector of eleven-dimensional super- 
gravity. We find new cosmological solutions that can dynamically explain the existence of 
three large spatial dimensions characterised by an universal asymptotic scaling behaviour 
and a large number of initially unwrapped dimensions. This type of solutions enlarge the 
possible initial conditions of the Universe in the Hagedorn phase and consequently can 
potentially increase the probability of dynamical decompactification from anisotropically 
wrapped backgrounds. 
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1. Introduction 

Understanding why we live in 3+1 dimensions is an intriguing and puzzling problem in 
theoretical physics. In the standard theory of general relativity the dimensionality of the 
Universe is an assumption and it cannot be derived dynamically or from a fundamental 
law. A suitable theoretical framework to deal with this question is provided by string 
theory which predicts the existence of extra dimensions. The general lore is to assume 
that these dimensions are effectively very small and consequently unobservable today. It 
is, however, fairly plausible to belive that this simple picture changes in the far past. As 
the Universe shrinks and gets smaller all directions should behave in the same footing and 
scale with similar sizes. Until definite physical evidence of the contrary is put forward one 
can consider the very early Universe to be higher dimensional and try to unveil a dynamical 
process responsible for the late asymmetry between large and small spatial directions. 

An important ingredient of string theory is the presence of extended objects. The 
idea that these objects can play a fundamental role in explaining the current number of 
dimensions was firstly suggested by Brandenberger and Vafa Q (see also |2|, ^, ^). In 
this proposal the background spacetime is considered to have (d+l)-dimensions with a 
spatial toroidal section and the dynamics is assumed to be driven by a gas of fundamental 
strings. Basically, all spatial directions are initially small, namely with a typical size 
of the order of the string length. As the Universe expands the energy of the winding 
modes that appear in the spectrum increases leading to a confining potential that stops 
the expansion. Then, to have a large Universe one has to devise a way to get rid of these 
modes. An efficient annihilation mechanism of string winding modes can only take place if 
the dimensionality of the spacetime is smaller than 4. This can be very easily understood 
on topological grounds. If the dimensionality of the spacetime is higher the probability of 
interaction of one dimensional extended objects is basically negligible. The picture is that 
of a Universe starting with all directions oscillating independently around the string length 
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scale until incidentally three of them get slightly larger than the others. Then, the process 
of annihilation is triggered letting the Universe grow with the appropriate dimensions. 

Recently, it has been proved that this mechanism also works if the dynamics of D- 
branes is taking into account The key argument is again topological because the 
probability of interaction of D-branes depends basically on the dimensionality of their 
worldvolumes. The larger the dimensionality the most effective the decay process is. In this 
picture, the dimensionality of the background spacetime is successively lowered as branes 
with larger dimensions are annihilated. Since the last winding modes to decay are those of 
the one-dimensional D-branes, an expanding Universe with three large dimensions and a 
hierarchy of small dimensions can be explained. Some particular aspects of the cosmology of 
brane gases have been studied by different authors [@, |, |, |ll], 11, ||, 0, |l8|. 
An interesting question that arises is whether brane gases can naturally stabilise the small 
extra dimensions at the string scale. In fact, this has been tested and verified in different 
set ups [|^, |2^, |2^, 24|. However, as it has been pointed out in |25], it seems quite 

hard to stabilise simultaneously the volume of the internal space and the dilaton without 
invoking new physics. 

The problem of dilaton stabilisation also brings a conceptual drawback. The existence 
of the dilaton itself is a consequence of assuming a priori that one of the dimensions is 
already compactified on a circle. Then, to have a proper mechanism of dynamical decom- 
pactification one should have started from a more fundamental theory. This was the pur- 
pose of the framework proposed in [26| for studying the cosmology of brane gases. Here, the 
starting point is the purely gravitational sector of eleven-dimensional supplemented with 
a gas of massless supergravity particles and a gas of nonrelativistic M2-branes wrapping 
anisotropically the spatial directions. Dimensions which are unwrapped become large at 
late times because of the absence of negative pressures that could suppress their expan- 
sion. Then, this scenario can explain a proper hierarchy of dimensions if there are initially 
three unwrapped spatial directions. Although a hierarchy of dimensions appears rather 
naturally the small dimensions do not stabilise in the relevant cases. As noted in a 
decompactification mechanism with isotropic wrapping can be also obtained if the inter- 
sections between M2- and M5-branes, which represent string-like degrees of freedom, are 
included in the dynamics. 

How plausible anisotropically wrapped configurations are depends on the thermody- 
namics of the brane gas close to the Hagedorn temperature. The analysis of this phase 
reveals that anisotropic wrappings with a low number of unwrapped dimensions are only 
compatible with a large initial volume of the Universe |27|. Basically, this is a consequence 
of the fact that for small volumes the rate of annihilation of branes increases and a larger 
number of dimensions can get unwrapped more rapidly. However, to put bounds on the 
initial conditions of the Universe and see whether a number of unwrapped dimensions is 
preferred at the time in which branes and antibranes freeze out new arguments are needed. 
A speculative possibility is to impose that the very early Universe is consistent with the 
holographic principle. This principle requires that the entropy inside a spherical volume 



must be bounded by the surface area 29 1. Estimating the entropy density close to the 



Hagedorn phase one can check that the initial state of the Universe should be rather small 
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and, then, initial brane configurations with a smah number of unwrapped directions are 
strongly disfavoured |27]. Nevertheless, this conclusion was obtained assuming isotropy of 
the background spacetime and should be taking very cautiously. 

The main purpose of this work is to investigated how the dynamics of a gas of branes 
changes if the gauge sector of eleven-dimensional supergravity is also taking into account 
and see whether some of the previous problems can be solved. There are not many works 
analysing the importance of gauge fields in the context of brane gas cosmology. In |lTl|| , 
it has been noticed that fluxes can be responsible for the existence of a subhierarchy of 
small dimensions in scenarios with M2-M5 intersections, and in [15|, it has been proved 
that a gauge field can dominate over the confining potential of the winding modes at late 
times and be responsible for the expansion of three large directions in a ten-dimensional 
supergravity background. 

Naively, one would expect that fluxes with a non-negligible dynamical contribution to 
the dynamics at late times will spoil the decompactification mechanism with anisotropic 
wrappings. Generically, a cosmology driven by a 4- form gauge field strength in a spa- 
tially flat background has seven expanding and three contracting spatial dimensions pO|| . 
Although a hierarchy of dimensions is dynamically created, it does not predict the right 
number of large dimensions. We will argue that in fact the dynamics of fluxes can be 
easily accommodated in this decompactification mechanism. Furthermore, it is shown that 
the presence of Euxes allows a new class of cosmological solutions with a large number of 
unwrapped dimensions which can account for three large spatial dimensions at late times. 
Solutions with this type of brane configurations make less stringent the initial conditions of 
the Universe in the Hagedorn phase and consequently enlarge the probability of dynamical 
decompactification from anisotropically wrapped brane gases. 



2. Brane gas dynamics in eleven dimensions 



We consider the eleven-dimensional = 1 supergravity with the fermionic sector frozen 
out. The degrees of freedom of this theory consist of a graviton and a 3- form gauge field 
A[3]. The effective action is given by pT| , 
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where R is the scalar curvature, F[4] = dA^^j is the field strength of the gauge field, and 
the eleven-dimensional gravitational coupling constant ku is given in terms of the Planck 
length by ^nf^ = (27r)^/]*^. In what follows, we will use Planck units in which In = 1. 
The last term is the Chern-Simons contribution and arises as a direct consequence of 
supersymmetry. In this work we only consider gauge fields with vanishing Chern-Simons 
term. For this class of solutions the classical equations of motion are simply. 



V/,F^^"^ = , 



(2.2) 
(2.3) 
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where T^^ is the energy-momentum tensor associated with the gauge field, 
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(2.4) 



and T^^ stands for the energy-momentum tensor of any other matter component. Note also 
that because the field strength is an exact form the dynamics must also obey the Bianchi 
identity V[pFf^i^af3] = 0- Since we are interested to investigate under which conditions three 
spatial dimensions grow large it is natural to consider spacetimes which are homogeneous 
but spatially anisotropic, 
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We further assume that the spatial dimensions have the topology of a torus so that the 
spatial coordinates have a finite range < x, < 1 and the total spatial volume is, 
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The non-vanishing components of the Einstein tensor for this metric are, 
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The equation of motion for the field strength (|2.3|) can be straightforwardly solved using a 



Freund-Rubin ansatz 133] (some cosmological solutions have been investigated in [33, 34 



30, 35| ). In this type of solutions the antisymmetric field strength is considered to be 
nonzero only on a 3-|-l-dimensional submanifold, say, 

j^t^uaf^ = l^=F{t) , (2.9) 

where indices run from to 3, 54 is the determinant of the induced metric on the subman- 
ifold, e^^°^^ is the ordinary Levi-Civita density, and the function F{t) is given by 



F{t) = /e-^4(i) ...g-AioW ^ 



(2.10) 



with / a constant of integration. The corresponding energy-momentum tensor for the 
gauge field is diagonal and its individual components can be compactly expressed as. 



1 



F{t] 



where the ten-dimensional object e* is —1 for i = 1,2,3 and +1 for i 
energy-momentum tensor corresponds to a fluid with energy density. 
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and anisotropic pressures, 

P'a = ^'Pc ■ (2-13) 

Note that as long as the gauge field has a dominant contribution the spacetime is naturally 
separated into RxT^xT^. The presence of the gauge field also introduces a new length 
scale into the dynamics given by the inverse of the initial vacuum expectation value of the 
antisymmetric field strength, 



la- 



(F^.a/^F^'^"^) '^'-llr\ (2.14) 



where Iq gives the typical length scale for the initial size of the Universe. The gauge field 
will be dynamically relevant if scales with or, equivalently, if the integration constant 
of the gauge field strength is of the order f l^. 

Apart from the gauge field we still need to specify the rest of matter sources before 
trying to solve Einstein equations (^]^). We assume that an important component of 
supersymmetric matter is present in the early Universe. This source of relativistic matter 
can be represented by a gas of massless particles, with energy density Pg and pressure 
p^. For simplicity we take the gas to be a homogeneous and isotropic perfect fluid with a 
radiation equation of state Pg = p^/W. The corresponding energy-momentum tensor is, 

(r,)^ = diag(-p„p„...,p,). (2.15) 

Note that because this energy-momentum tensor is covariantly conserved, the energy den- 
sity of supergravity particles scales with the total size of the Universe as, 



Vo 



11/10 



Ps = Pi [y ) , (2.16) 

where p° and Vq are, respectively, the energy density and the spatial volume at some given 
time, to- 

The second source of matter in this model is a gas of M2-branes wrapped on the various 
cycles of the torus. This gas can be characterised by a matrix of wrapping numbers Nij, 
where elements with i < j represent the number of branes wrapped on the (ij) cycle while 
elements with i > j represent the number of antibranes in the same cycle. The elements 
of the diagonal are irrelevant and can be chosen equal to zero. Since we assume that 
the number of branes and antibranes is the same for each cycle the wrapping matrix is 
symmetric. For the later discussion it is useful to classify the spatial dimensions into three 
types. A direction i is said to be unwrapped if Nij = for all j. Fully wrapped directions 
have all Nij nonzero except those that correspond to an unwrapped direction. A direction 
for which some of the components Nij are zero for values of j corresponding to a not 
unwrapped direction is referred to as partially wrapped. The actual values of this wrapping 
numbers should be provided by a, still lacking, theory of thermal and quantum fluctuations 
in the very early Universe. For our purpose, and until our theoretical understanding of 
this underlying theory is improved, we will simply take these numbers as random integers. 
Note that although the theory also supports the existence of M5-branes, their dynamics is 
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not taking into account because they decay very rapidly in eleven dimensions and cannot 
produce any relevant effect at late times. 

A single M2-brane is describe by the Nambu-Goto action, 

Sm2 = -T2 [ d?a (2.17) 

JM-i 

where the coordinates a"'{a = 0, 1,2) parametrise the three-dimensional worldsheet mani- 
fold A^s spanned by the brane and 53 is the determinant of the pull-back onto ^A■i of the 
eleven-dimensional bulk metric g^y. The surface tension of the brane is a fixed parameter 
given in Planck units by |36|, 

Following [^] we will ignore excitations on the brane worldvolumes and assume that the 
branes are non-relativistic. Under these conditions, the non-vanishing components of the 
stress tensor uniformly averaged over transverse directions for a brane gas with wrapping 
matrix Nij will be given by, 

Pb = -(Tjt = J^:;^E^''e''N,,, (2.19) 



k^l 

1 



p; = {Tj, = -^_A_^e^'e^'= (TV^, + iV.fc) . (2.20) 



where we have introduced the total energy density, pg, and the anisotropic pressures, p^, 
respectively. The above M2-brane action assumes that the brane is moving in the eleven- 
dimensional background without interacting with the gauge field. Nevertheless, it is easy 
to check that under our assumptions coupling terms of the form 



T2 / ^[3] , (2.21) 

where A[^j{x{a)) represents the pull-back of the eleven-dimensional gauge field 3- form onto 
A43, do not contribute to the dynamics and can be safely neglected. 

Now, we have to insert all the matter sources of energy-momentum, T^'^ ,T^^ ^ and 
T^'', into the right hand side of Einstein equations (^.2|). After some algebra the time 
component of the equations can be expressed as, 

10 

Y,\kXi = 2Kl^{p,+Pg+p^) , (2.22) 

k^l 



and the spatial components as, 

^ ■ ^ ■ i = l,---,10, (2.23) 



Here, H is the mean Hubble parameter which represents the rate change of the total spatial 
volume of the Universe and is given in terms of the metric components by, 

10 

H = -Y.\^. (2.24) 
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One can readily observe that now there are two possible source of anisotropy in the cos- 
mological evolution. The one coming from the brane gas and the one from the gauge 
field e:^ Pq ■ A quantitative way of measuring the anisotropy of a spacetime is by means of 
the shear scalar which can be defined for our metric ansatz as, 



1 ^° 2 
^'(*) = qE(^^-^) • (2-25) 



10 

9 

Obviously, this object is zero only if all the expansion rates are equal. That the shear 
scalar can be nonzero in our model can be explicitly seen by comparing the evolution of 
two different expansion rates, 

A, - A, + im{K - A,) = K?, [{e^ - e^)p^ + (p^ - p>^)] . (2.26) 

From this expression one easily observes that both the gauge field and the gas of M2-branes 
are the sources of a nonisotropic evolution of the Universe. The pressures, p^, exerted by 
an anisotropically wrapped brane gas are always nonpositive and, then, their effect is to 
suppress the growth of the scale factor in the dimensions with nonzero wrapping numbers. 
An asymmetric wrapping, in the sense that some spatial dimensions are unwrapped and 
then differential pressures are allowed, will lead to an anisotropic expansion. On the other 
hand, the effect of the gauge field pressures, = eVc depends on the sign of e*. For those 
with e* = +1, the size of the spatial dimension i is enhanced and for those with e* = — 1 
suppressed. Two spatial dimensions i and j will have different expansion rates if i= ■ 
Hence, it is of fundamental importance for understanding the asymptotic cosmological 
evolution to see which energy components dominate at late times for different initial brane 
configurations. This is better done by introducing a new set of dimensionless functions of 
cosmic time, 

^.W = ^fe> ^oi.t) = ^%- (2.28) 



Here, , Vtg , are the ordinary density parameters for each matter component (super- 
gravity gas, brane gas, and gauge field, respectively), and ilo- the fractional contribution 
of the shear to the expansion of the Universe. With these definitions the constraint ( 2.22| ) 
becomes simply, 

Q.„ + VLg+Vtg + n^^ = l. (2.29) 



3. Dynamical features at late times 

Let us start the analysis of the late time dynamics by describing some general properties. 
Without a gas of branes and no gauge degrees of freedom, flat and Kasner spacetimes are 
exact vacuum solutions with = 0. When the energy density of the gas of supergravity 
particles is nonzero, the cosmological dynamics is that of an isotropic eleven-dimensional 
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radiation dominated universe with a growing scale factor ~ t^/^^. On the other hand, 
for pg = and a gas of branes wrapping all the spatial dimensions isotropically an exact 
solution exits with e'^ ~ t^^^ [^] . It is easy to show that in general the energy density of the 
supergravity particles cannot fully dominate at late times if a brane gas with a nonisotropic 
wrapping matrix is present. As a consequence, the final fate of the Universe in this case is 
always anisotropic and a hierarchy of dimensions can always be created. Finally, if 7^ 
and all other energy components vanish, one can find an asymptotic power-law solution, 

(3.1) 

This clearly means that the dynamical effect of the form field is to drive a cosmological 
evolution of the Universe with 3 contracting and 7 expanding spatial dimensions. Then, if 
the asymptotic evolution of these type of cosmologies is fully dominated by the gauge sector, 
even in the presence of wrapping branes, the correct number of large spatial dimensions 
cannot be explained and the present decompactification mechanism would be completely 
invalidated. Fortunately, as we will see shortly, this is not the case. Another important 
point to stress from this scaling analysis is that a higher-dimensional Universe dominated 
by a 4-form field has subvolumes that shrink. Consequently, unless quantum effects prevent 
the contracting dimensions from reaching a zero size, a physical singularity appears in the 
future. Furthermore, one can find exact solutions in which this singularity occurs at a 



finite proper time ||30[]. Note that, this picture can qualitatively change if the underlying 



spacetime has a submanifold with spatial curvature. 
3.1 Cosmological evolution without fluxes 

First consider the dynamics without fluxes. It is important to remark that all the comments 
we are going to make in this section are also applicable to those situations in which the 
fluxes are negligible at late times but not necessarily at intermediate stages. 

We take initially a Universe filled with a gas of supergravity particles and a gas of M2- 
branes described by a wrapping matrix with mi unwrapped, m2 partially wrapped, and 
m-3 = 10 — mi — m2 fully wrapped spatial dimensions. Generically, we will refer to such 
brane configurations as a brane gas with mi-m2-m^ wrapping. In general, the solutions of 
the Einstein equations (|2.23| ) present an universal power-law behaviour at late times. That 
is, the sizes of all the spatial dimensions grow with a power law which depends solely on 
the wrapping matrix and not on the initial conditions [p6| . These attractor solutions have 
a generic analytical form given by. 



t'^ for i = 1, • • • , mi (unwrapped) 

for i = mi + 1, • • • , mi + m2 (partially wrapped) (3-2) 
for i = mi + m2 + 1, • • • , mi + m2 + m3 (fully wrapped) 



As discussed in [26|, a physically relevant hierarchy of dimensions is produced only if 



mi = 3. For instance, a brane gas with a wrapping matrix of type 3-3-4 has, 

a = ^, P = i = Yr ^^-^^ 
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and another with wrapping matrix of type 3-4-3, 



These two examples illustrate that a gas of branes with the appropriate wrapping can 
support naturally a faster growth of three spatial dimensions because a is always greater 
than P and 7. This is because the brane gas exerts no pressure on unwrapped directions and, 
consequently, cannot suppress their expansion. For partially and fully wrapped dimensions 
the pressures are negative and cause the opposite effect. One also observes, that there 
exist configurations in which /3 > 7 and a sub-hierarchy in the sizes of small directions can 
also be formed. This feature occurs in configurations with wrapping parameters satisfying 
m2 > m^. In fact, one can further show that the powers only depend on the total number 
of spatial dimensions (in our case 10) and the wrapping matrix parameters mi, 1712, and 
m3 (the explicit expressions can be found in p6| ). 

An account of the full dynamics can be obtained by solving Einstein equations nu- 
merically (see Fig |l]). Since we do not expect supergravity to hold at very early times, 
the strategy we have followed is to evolve the system from the Planck time onwards given 
initial random values for Aj and Aj. For the nonzero components of the wrapping matrix 
Nij we take random integers chosen from the set {1, 2, 3}. We have checked for consistency 
that the final expansion rates for each spatial dimension are independent of these choices. 
The constraint ( 2.22| ), or ( 2.29| ), is only used to compute the initial value of the energy 



density of the supergravity gas, p°, and to check the numerical accuracy of the output. In 
our particular numerical implementation, the universal powers found in all the examples 
agree with those obtained analytically with less than 1% error. In general, the dynamics 
presents three stages. A short initial phase in which any initial anisotropic expansion is 
rapidly diluted, a much longer intermediate phase of isotropic expansion, and a final phase 
in which the unwrapped dimensions grow with an expansion rate larger than that of the 
rest of spatial dimensions. As one can observe from the plots depicting the contributions to 
the constraint ( 2.2S| ), the shear, which represents the degree of anisotropy of the Universe, 



decreases in the first phase, is negligible in the second, and grows again until it scales with 
the mean expansion in the third. The energy density of the supergravity particles domi- 
nates the expansion of the Universe in the first two stages and the energy density of the 
brane gas in the final stage. The time at which both contributions are of the same mag- 
nitude marks the transition between the last two phases. For the simulations represented 
in Fig. H this time is approximately t ~ 10^ but the actual number depends on the initial 
conditions and the wrapping matrix chosen. 

In conclusion, the cosmology of brane gases in a (low-energy) M theory context opens 
the possibility of explaining the number of spatial dimensions observed today for those con- 
figurations characterised by a wrapping parameter mi = 3. In the following we see that this 
conclusion is still true even if fluxes have a significant contribution at late times. Further- 
more, we present new cosmological solutions with a larger number of initially unwrapped 
directions which support asymptotically a large four dimensional spacetime. 
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Figure 1: Cosmological evolution without fluxes. Graphs on the left are for a brane gas with a 
wrapping matrix of the type 3-3-4 and those on the right of the type 3-4-3. Plots on top represent 
the time evolution of the size (scale factor normalised by 27r) of all the spatial dimensions: solid 
curves are for unwrapped dimensions, dashed curves for partially wrapped dimensions, and dotted 
curves for fully wrapped dimensions. Plots on the bottom depicts all the contributions to the 
expansion of the Universe as a function of cosmic time. The solid line represents , the dashed 
line fig , the dotted line fig- , and the dotted-dashed line il^ . The thick solid line is the sum of all 
the contributions and serves to check the accuracy of the numerical computation (2.29). 



3.2 Cosmological evolution with fluxes 



Now we are interested to investigate the influence of fluxes on the late-time cosmological 
dynamics of a Universe fllled with a gas of M2-branes. The flrst difficulty one faces when 
fluxes are present is that the wrapping matrix describing the branes gas and the gauge field 
induce different splittings of the spacetime. As we have seen, the type of solutions for the 
4- form field strength we are considering naturally separates the spacetime into RxT^xT^ 
and the gas of M2-branes into RxT*"! xT™^ xT™'^. In spite of that, the universal power- 
law scaling behaviour at late times we have described in the previous section, although 
slightly modified, is not lost. 

Given a brane gas configuration with wrapping matrix of type mi-m2-m-i the late 
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cosmological evolution can be cast into the form, 



(3.5) 





for i 




for i 




for i 




for i 



I,-- - ,3 (unwrapped) 
4, • • • , mi (unwrapped) 
mi + 1, • • • , mi + m2 (partially wrapped) 
mi + m2 + 1, • • • , mi + m2 + m^ (fully wrapped) 

where we have chosen the first three directions to be those picked by the Preund-Rubin 
ansatz for the gauge field. As seen in Fig. for configurations with a wrapping matrix 
3-m2-m3, the energy density of the gauge field and the brane gas both have a significant 
contribution at late times. In this case (note that there are no directions with a_|_) the 
scaling behaviour is given analytically by, 

a- = ^, /3 = 7 = ^, (3.6) 

when m2 < m^, and, 

^ _ 5m2m3 - 14 ^ _ 5m3(m2 - 3) + 7 ^ _ 5m2(m3 - 3) + 7 
5?7i2m3 + 49 ' 5m2ms + 49 ' 5771,2^3 + 49 ' 

when m2 > m^. As a particular example of the last case, one has 

a_ = — , /?=— , 7=—, (3.8) 
109 ' 109 ' ' 109 ' ^ ^ 



for m2 = 4 and 7713 = 3. Comparing these results with ( |3.3| ) and ( |3.4D one readily observes 
that, although in both types of configurations the expansion rates of the unwrapped direc- 
tions are decreased and those of the partially and fully wrapped are increased, the inEuence 
of Euxes at late times does not destroy the hierarchy between large and small dimensions 
(a_ is still greater than /5 and 7 in all the cases). Note as well that even though the 
energy density of the gauge field fully dominates at early times, the energy density of the 
gas of branes very quickly takes over and the collapse of three dimensions is avoided. The 
appearance of a physical singularity in this early phase of the evolution is then prevented 
at the classical level. 

A brane configuration of particular interest is that of a brane gas with 6 unwrapped 
dimensions; that is, configurations with wrapping matrix of type 6-m2-ms- Without fluxes, 
there are only four possible configurations leading to an anisotropic evolution. For all the 
configurations with m2 < m^ one has the analytic solution ( ^.2| ) with, 

4^1 . ^ 

« = YY' ^ = ^ = -n- ^^-^^ 

On the other hand, for the configuration with wrapping 6-3-1, which is the only one with 
m2 > m^ supporting a nonisotropic expansion, one obtains the scaling, 

« = |, P = Yr ^ = -1- (3.10) 

Obviously, if the dynamical effects of the gauge field are not taking into account these 
configurations, although providing a hierarchy among different dimensions, do not predict 
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Figure 2: Cosmological dynamics of a brane gas of wrapping type 3-3-4 (left) and type 3-4-3 
(right) with fluxes. The identification of plots and lines is the same as in Fig. |[ 

the right number of large spatial directions. Finally, it is important to note that, contrary 
to what happens when mi = 3, the energy density of the supergravity particles is always 
not negligible at late times. 

Let us now consider the dynamical effects of the gauge degrees of freedom in config- 
urations with mi = 6. Two numerical examples are compared with the no flux case in 
Figs. 1^ and ^. An interesting point is that all the energy components (fi^, ilg, fi^) have 
a significant contribution to the cosmological expansion at late times. As one can also 
observe, these solutions represent a new family of configurations that permits the growth 
of three spatial dimensions. Although for the configuration with wrapping matrix 6-3-1 the 
hierarchy of the larger dimensions is not quite obvious from the plot obtained numerically 
(Fig. ^) , it can be seen analytically that the asymptotic behaviour of the solution has the 
form (|3.5| ) with, 

"- = 44' ^=44' 

"+ = 44' ^ = "44- ^^-^^^ 
The values obtained in our numerical computation agree with the above analytical result 
up to four decimal places. Comparing with the case without fluxes, /3 is larger and 7 
smaller. That means that the three partially wrapped dimensions grow faster and the fully 
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Figure 3: Cosmological dynamics of a brane gas of wrapping type 6-1-3 without (left) and with 
(right) fluxes. The identification of plots and lines is the same as in Fig. |l]. 

wrapped dimension contracts also faster. On the other hand, the six unwrapped dimensions 
are separated into two groups, of three dimensions each, with independent expansion rates 
defined by a_ and a+, respectively. Since a_ < a+ the appropriate hierarchy for explaining 
the number of spatial dimensions is produced. On the other hand, for m2 < ms one finds 
analogously the asymptotic power-laws (see Fig. |3|), 

3 5 1 

Again a+ > a_ > /5, 7 and a hierarchy with three dimensions growing larger is obtained. 
Consequently, brane configurations with nii = 6 unwrapped dimensions can naturally 
explain why three large spatial dimensions are decompactihed at late times if the gauge 
sector of eleven-dimensional supergravity is turned on. This new family of solutions enlarge 
the possible wrapping configurations at the end of brane-antibrane annihilation in the 
Hagedorn phase. It is important to emphasis that the interplay between both the brane 
gas and the flux dynamics plays a fundamental role in getting the correct number of large 
dimensions. 

4. Conclusions 

We have illustrated how the cosmological evolution at late time of a gas of M2-branes 
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Figure 4: Cosmological dynamics of a brane gas of wrapping type 6-3-1 without (left) and with 
(right) fluxes. The identification of plots and lines is the same as in Fig. |l]. 



within a low-energy limit of M theory is modified when the gauge fields of the bosonic 
sector are taking into account. 

We have seen that fluxes respect the hierarchies among different spatial dimensions 
introduced by an anisotropically wrapped brane gases at late times. We have also found 
new solutions that can explain the actual number of spatial dimensions of the Universe 
which are characterised by a brane gas configuration with a large number of unwrapped 
dimensions. These solutions appear as far as the gauge field strength is sufficiently strong 
to have a significant contribution to the total cosmological expansion at late times. On 
thermodynamical grounds one should expect that these brane configurations can be origi- 
nated at the end of the Hagedorn phase from smaller initial spatial volumes of the Universe 
than configurations with lower unwrapping numbers. Including gauge fields into the dy- 
namics increases the possible initial conditions for the Universe and then the probability 
of obtaining decompactification from anisotropically wrapped spacetimes. This makes less 
severe the fine-tuning problem posed by this mechanism when fluxes are not considered. 

A difficult issue is to assess whether the constraints on the initial size of the Universe 
imposed by the holographic principle are significatively modified when gauge degrees of 
freedom are included. For the type of solutions we have studied the assumption of isotropy 
as considered in [27 \ is certainly not acceptable and a fully anisotropic analysis of entropy 
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bounds is absolutely compulsory. In principle, the presence of a new physical parameter 
representing the strength of the gauge field opens the possibility of getting less restrictive 
conditions on the physical volume of the Universe in the Hagedorn phase. 

Unfortunately, as in the case without fluxes, the internal dimensions do not stabilise 
for the physically interesting cases. In general, it seems rather difficult to get stabilisation 
within this mechanism without introducing additional physics. 

Finally, it would be certainly interesting to investigate if there exist solutions with an 
even larger number of unwrapped dimensions that could as well explain the number of 
spacetime dimensions. A full classification of solutions and configurations will determine 
how generic is this anisotropic mechanism of decompactification. This analysis will further 
reveal the connection with the mechanism suggested in jl^ within the context of ten- 
dimensional type IIA supergravity. This lower dimensional theory arises as a result of 
the compactification of eleven-dimensional supergravity and, thus, in principle, one should 
expect a close relationship between the cosmological solutions in both frameworks. 
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